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1. INTRODUCTION 


1 . 1 WHENCE SWIRL 

Swirling of a liquid in a container may, in some applications, 
prove to be a better way of managing liquids in space than 
rotating the entire container. A transient study by the National 
Bureau of Standards in Gaithersburg Md. I 1 3 showed that swirl 
could be started quickly by the injection of a relatively high 
velocity jet of liquid tangentially into the body of liquid. One 
unknown still remained: the quantity and velocity of jet flow, or 

mechanical power of a solid impeller required to maintain a given 
radial acceleration ("G" force). 

The crux of this problem is to determine the rate of rotational 
energy dissipation by wall friction. 

1 . 2 RECAP 

The part of this study which ended in December ' 86 gave a first 
approximation of the jet requirements to maintain a given "G" 
force in a tank. A sample calculation was made for a hypothetical 
5,000 liter liquid helium tank with a radial acceleration of 
O.OIg. The resulting estimation of jet flow was 3(10)~ 5 m 3 / s. If 
the NBS (Boulder) centrifugal pump were used for transferring 
liquid helium from this tank, the jet flow just stated would 
amount to about 22% of the pump output. Since the required jet 
flow is a significant portion of the total pump output, the 
accuracy of calculation could affect a decision as to the 
feasibility of the jet method, or the suitability of a particular 
pump. 

The present study investigates the possible additional effects of 
axial variation of tangential velocity and secondary (radial and 
axial) flow components. 

1 . 3 SUMMARY OF THE PRESENT STUDY 

1.3,1 Part 1 (under the assumption of negligible secondary flow) 

Since the rotational energy dissipates along the cylinder walls 
and end closures, it seemed likely that the tangential velocity 
and angular momentum would diminish appreciably with axial 
distance from the input. Therefore, the first investigation 
concerned the effect of axial variation of angular momentum 


-I - 



assuming negligible radial and axial velocity components. This 
assumption together with a simplified turbulent viscosity reduces 
the tangential equation of motion to a two dimensional ( radial 
and axial) linear equation solvable in closed form. The solution 
produces both the tangential velocity radial profile and an axial 
decay factor. 

The linear theory predicted that the velocity would decay to 1/10 
of its input plane velocity in an axial distance of one radius. 
Simple visual experiments with a spinning horizontal disk in a 
vertical cylinder revealed a much different behavior 
characterized by: 

1. A thin boundary layer on the spinning disk, cylinder walls, 
and bottom plate. 

2. Very little axial decay of tangential velocity in the body 
of liquid. 

3. The fluid body spinning at 1/2 to 1/6 the rotational 
speed of the disk (depending on the disk-bottom spacing). 

4. Noticeable outward radial flow near the spinning disk, 
downward flow along the cylinder wall, inward flow at 
the bottom, and upward flow at the center. 

To test the possibility that the linearization was at fault, the 
last attempt at "pure rotational flow” employed Prandtl* s mixing 
length hypothesis for turbulent shear and required a numerical 
solution of the equation of motion. The resulting axial decay 
was only slightly less severe. 

1.3.2 Integral momentum approach 

Even though the radial and axial velocity appear insignificant, 
these secondary flows have a large effect as evidenced by the 
failure of the "pure rotation" theory to describe the actual 
flow. Therefore, the second approach applies a published method 
developed to predict the frictional torque due to the housings of 
rotating turbine or pump impellers (a three dimensional flow 
problem). For the swirl problem it was necessary to extend the 
theory to account for the cylinder walls and the existence of a 
free interior surface in the absence of gravity. The 
relationship between the cylinder length, drag coefficient, and 
the difference in speed between a spinning disk and the liquid 
was demonstrated in a general way by visual experiments. 

The sample calculation for the 5,000 liter liquid helium dewar, 
repeated using the new drag coefficients, resulted in a much 
larger jet flow requirement. In terms of the NBS centrifugal 
pump, about 55% of its flow would be required to maintain 0. Olg 
radial acceleration instead of the 22% calculated earlier. 
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2.0 ANALYSIS OF SWIRLING FLOW 


2. 1 AXIAL VARIATION OF TANGENTIAL VELOCITY ASSUMI NG NEGLI GI BL E 
SECONDARY FLOW 


In the first part of this study it is assumed that radial and 
axial velocity components are negligible, ie. the flow is 
assumed to be purely rotational, but with tangential velocity 
varying in both radial and axial directions. 

It should be noted that the final results cast doubt on the "pure 
rotation” hypothesis, and the work done under that hypothesis is 
only briefly outlined. Details are given in Appendix A. 

2.1.1 A Linear Solution 

The assumptions are: 


1. Steady state 3/3t = 0 

2. Pure rotation v P =0, v x =0 

3. Axial symmetry 3/39 = 0 

4. No external forces g«=0 


By the simplification of using an average value of the turbulence 
factor, I 3m/ 3z I / I 3m/ 3r I = S 2 , rather than considering it a 
variable, the 9 equation of motion reduces to a linear equation 
in r and z, solvable in closed form. The z equation is 
exponent i al: 

m( r, z) = m( r, 0) exp( z/Zd ) , where (1) 

2d is the z decay constant whose numerical value comes 
from the solution of the equation in r. 

The equation in r is a Bessel equation with the following general 
solution: 

. [c,j, (**) ♦ C,r, (**>1 (2) 

run L v Zd^J 

Ji is the Bessel function of the first kind, first order. 

Yi is the Bessel function of the second kind, first order. 

The values of the constants Ci, C 2 , and Zd are dependent on the 
fullness of the cylinder (ri/r o =0 being full and ri/r 0 -1 being 
empty), and the Reynolds number. These constants are evaluated 
in tables Al and A2. The resulting values of Zd/r 0 being near 
0.4 means that, in an axial distance of one radius, the 
tangential velocity would decay to less than 1/10 its value at 
z=0. A very simple visual experiment shows clearly that this is 
not true; in fact, it appears that Zd/r 0 should be at least 2. 
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The calculated variation of Z from equation (2) reveals a spike 
near » even though a constant value was used in the derivation. 
To test the possibility that the linearization (S=const.) causes 
the unrealistically small Zd/r 0 , the non-linear turbulent equation 
was solved numerically. The resulting z<j/r 0 of 0.46 is hardly 
an improvement, and so it appears that the linearization was not 
the problem. 

2. 2 ADAPTATION OF A MOMENTUM INTEGRAL METHOD FOR TURBULENT DRAG 
DUE TO SWIRLING FLOW IN A HOUSING 

The unrealistic decay is apparently caused by neglecting the 
secondary circulation which transfers momentum by large scale 
movement of fluid containing the momentum. This circulation, 
which takes place mostly within the boundary layers (or Ekman 
layers [ 4] ) , can be very powerful compared to ( small scale) 
turbulent momentum transfer. Greenspan [41 states: 

secondary flow can significantly alter the primary motion 

through a slight redistribution of angular momentum and 
vort i ci ty" . 

The drag on a pump impeller or turbine wheel in a housing is a 
problem similar to the present one in which secondary flows play 
an important role. Schul tz-Grunow employed the momentum integral 
method as described in Greenspan C 43 and Schlichting (33 and 
arrived at an expression for drag coefficient on a spinning 
disk in a housing. The connection with the present problem is 
that the drag on the rotating liquid, which we require, is the 
exact balance of the drag on the spinning disk. The spinning 
disk may either be a real or an imaginary device. 

The adaptation of the Schul tz-Grunow theory to account for an 
elongated cylinder and the possibility of free surfaces is 
detailed in Appendix B . The result is the following expression 
for frictional torque on a spinning disk (or torque of the liquid 
on the housing). Equations ( B5) with ( B9) : 

To = 0. 041 ( 2/ f ) 1 ' s puo 2 r„ 5 ( 1 / Uo ) 7 ' 5 (l -( ri /r. ) 2 3 ' 5 ) Re 0 * 1 ' 5 

( 3) 

Reo=u)or 0 2 /v/ 

The remainder of the problem consists of determining the ratio 
Ui/U o, and the factor f, both of which are functions of the 
geometry . 

Appendix B contains the details of the modifications of the 
Schul tz-Grunow equation to include: 

1. An elongated cylindrical section 

2. A possible free interior surface 

3. A possible solid interior concentric cylinder wall 
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(as in proposed ground experiments) 

4. A possible free top surface as in normal gravity 
swirl in a partially filled vertical cylinder. 


The rotational frequency ratio Uf/Uo is worked out for two cases 
which cover the four possibilities listed above: 

Case 1 . For a partially or completely filled cylinder in zero 
gravity having a free surface (ie. no torque) at ri , from (B15) 
and ( B1 6) : 


Ud /Uf = (F( L) +l] w 7 + 1 


where, 


F ( L) 


23 L r 

5 r 0 


r o c 
r o 


■] 


i a / 3 


1 - ( r i / r o ) 


2 3/5 


( 5) 


Case 2 . Swirl in an open top vessel in a normal gravity 
environment. There will be no torque on the top end (Te = 0). 

The following also allows for a possible inner fixed cylinder as 
in a proposed swirl test apparatus: 


Ud / Uf 


4/7 

[f( L) (l +( r i / r o ) 1 a 7 5 JJ + 


( 6 ) 


The factor f arises from the increase in the boundary layer 
thickness due to the added frictional surfaces. The quantity 
(2/f) l/5 can be considered unity except for very long cylinders. 


f = 2 + fc. 


( 7) 


where fc is the cylinder contribution. Equation ( B1 9) : 

r = 1 8 r Uf / Up 1 7/4 L r r o c i 15/5 1 

c " 5 Ll-Uf/UoJ r 0 r 0 1 -( r t / r« ) 1 a / 5 


Appendix B also presents experimental data which gives some 
credence to the modification. 

Comparison With Previous Drag Calculations 

Returning to the example of final report 12/86: 


5,000 K cylindrical tank 

Arbitrarily specified radial acceleration O.OIg 


Total length = 2L 

ri 

r a 

Uf = v 8 / r 


2. 94 m 
0. 635 m 
0. 735 m 
0. 393 s" 1 


r o c = r o 
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ri/r*=0. 864 
L/r, =2 

Fluid: normal liquid helium 

v=w/ p = 8. 28(10)"® m 2 /s 
p = 1 45 kg/ m 3 

Re 0 =1 . 66( 1 0) ® 

F(L)=18.8 (equation (5)) 

Ui/Wo=0. 154 (equation (4)) 
f = 2 . 79 ( equation ( 7) ) 

T = 2T 0 =0. 1 28 kg m 2 /s 2 

If an actual spinning disk were used to generate the swirl: 

Power to the spinning disk = T(Jo-0. 332 w. 

The frictional torque calculated in the 12/86 report was 0.0488; 
the new value is 2. 7 times the old for the same conditions. The 
jet flow requirement, being approximately proportional to the 
frcitional torque will also be 2.7 times the previous value, or 
57% of the NBS pump output to maintain O.OIg radial acceleration 

The torque also varies with the degree of fulness of the vessel 
as follows: 



r i / r o 

Ud/(iI( 

T 

Full 

0 

4. 77 

0. 1 39 


0. 2 

4. 77 

0. 1 39 


0. 4 

4. 80 

0. 1 39 


0. 6 

4. 97 

0. 1 37 

Example 

0. 864 

6. 51 

0. 1 30 

Empty 

1 . 0 

CO 

0 


3. 0 CONCLUSIONS 

The pure rotation assumption leads to a predicted sharp decrease 
in tangential velocity with axial distance from the input plane. 

A simple experiment with swirling flow in a cylinder demonstrates 
that axial decay of velocity is actually very much smaller. 

The unrealistic decay produced by the pure rotation assumption 
indicates that radial and axial circulation within the boundary 
layers have a strong effect on the momentum distribution and on 
the frictional drag. The published theory of Shul tz-Gr unow 
treats a similar subject of turbulent rotating flow with radial 
and axial circulation in a housing. This theory, modified to 
match the present application, predicts nearly three times the 
drag calculated previously for a sample case. A far thinner 
boundary layer is the reason for the increased drag. Simple 
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visual experiments lend credence to the modified theory. 


A swirl inducing jet flow was also three times as large as 
originally thought (57* the NBS pump’s output for the example 
5,000 R tank) and this would mean: 

1. A larger pump than the NBS pump would be required for 
tanks the size of the example, or 

2. A smaller radial acceleration than O.OIg would be 
achi e ved. 

Since frictional torque is proportional to the vessel surface 
area, the NBS pump suffices for a smaller tank. 

Other means of producing swirl should be considered. The sample 
calculation shows that the mechanical power to achieve O.OIg in a 
5,000 8 tank is only 0.34 watts. 


4. 0 WHITHER SWIRL 

The Schultz-Grunow theory underpredicts his own experimental data 
by 17*, and the crude experimental observations reported here 
indicate that the modification for long cylinders may not be 
exactly correct. Therefore, it seems that theorizing has gone 
about as far as it should go whithout accurate experimental 
measurements of rotation rates and drag due to turbulent swirl 
in a scaled configuration similar to the application. Normal 
gravity experiments with a spinning disk, both at room 
temperature and with liquid helium, would be beneficial. 

A good check of jet flow rate to produce a given swirl velocity 
could be performed at room temperature with an apparatus only 
slightly more sophisticated than the one used to visualize the 
rotational velocity. Experiments of this type would be a useful 
supplement to the planned helium swirl tests. 

The findings of this swirl study regarding secondary flow, as 
well as the extensive studies of Greenspan, have a bearing on the 
spin-up of liquid in an impulsively started cylinder. It seems 
that spin-up time could be greatly reduced by intentionally 
stimulating radial and axial flow during spin-up. Design and 
testing of suitable vanes to accomplish this would seem 
desirable. 


5 NOTATION 


C coefficient in Prandtl' s mixing length hypothesis 

Ci, C 2 , C 3 , C 4 constants of integration 



Cf 

c, 

f 

F(L) 

F» 

g 

g. 

J i , Jo, Yt, Yo 

R 

L 


m 

m* maximum 

p 

r ) r i ) r« , r o j 

Re, Ref, Reo -- 

t 

T, T c r i , To, Te 


V, Vr , V, , V,, 

Z , Z D 


coefficient of the expression for drag coefficient 

constant defined in equation ( B5) 

drag ratio defined in equation ( B2) 

function of geometry defined in equation ( B16) 

wall shear function defined following table A1 

acceleration of gravity 9.8 m/s 2 

external force per unit mass 

- Bessel functions used in equations ( A7) and ( A1 3) 

Prandt 1 ' s mixing length, m 

half length of a cylindrical vessel, m 

angular momentum per unit mass, kg m/ s 

angular momentum per unit mass in momentum profile 

pressure, Pa 

r oc radius, internal, at maximum momentum, 

outer, outer-cylinder, m 

Reynolds number, based on fluid velocity, 

based on disk velocity 

time, s 

Xj torque, on cylinder, on disk, 

on end disk, on internal cylinder, kg m 2 /s 2 

v, . velocity, radial, axial, tangential, 

tangential at maximum momentum, m/ s 
axial coordinate, axial decay length, m 


5 boundary layer thickness, m 

q tangential angular coordinate 

|j T viscosity, turbulent viscosity, kg/m-s 

kinematic viscosity = u/p, m 2 /s 

£ ratio of momentum gradients, equation ( A5) 

p density, kg/m 3 

T, Tr i , Tii — shear stress, r-9 component, 0-2 component, kg/m-s 2 
u, Ud , Uf rotational frequency, of disk, of fluid, s' 1 
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appendix a 


AXIAL VARIATION OF TANGENTIAL VELOCITY 
ASSUMING NEGLIGIBLE SECONDARY FLOW 

A LINEAR SOLUTION FOR PURE RO TATI ON . IN TWO DIMENSION S 


A. 1 

The 9 component equation of motion in cylindrical c 
in ter ms of shear stress [11: 


oordi nates is, 


'3vi 9vi + X,. Iv,. + ♦ Pg. 

_3t 3r r 39 r 3z J 30 

Pi d f 2_ •) . L iki + 

' Lr 2 8r r 39 3z J 

The following assumptions will be made: 

1. Steady state . 3/9t ”° 

2. Pure rotation Vr= °a/ooIn 

3. Axial symmetry 

4. No external forces g * " 


( A1 ) 


Then ( 1 ) becomes: 


K f (r 2 Tr • 3 ♦ 

r 8r v 


8t 


8z 


= 0 


These equations apply to turbul 
when appropriate expressions ar 
The expressions for turbulent s 
complex, even differential equa 
Spaulding [ 21 point out that Pr 
( MLH) is still useful at least 
is relatively simple it is used 
development of MLH for flow ove 
angular momentum transport for 
substitutes torque for force, t 
rotating flow result: ( same ass 


ent flow as well as 1 
e used for the shear 
hear stresses can bee 
tions in themselves, 
andtl' s mixing length 
for approximations, a 
here. If, in Prandtl 
r a flat plate, one s 
linear momentum trans 
he following expressi 
umptions as above) 


( A2) 

aminar flow 
stresses, 
ome extremely 
Launder and 
hypothesi s 
nd since it 
• s intuitive 
ubstitutes 
port and 
ons for 


• [ u - “’■] r 


8m 
8z ’ 


C pX 2 8m 
= r 8z: 


Ut 0 


dm 

8r 


where 


( A3) 


[ u " UT# 1 r 8? + 2Um/i 

C£il| 


where 


( A4) 


In (A3) and ( A4) X is the mixing length, generally taken to be 
the distance from the nearest surface, and c is an empirical 
constant In what follows X and c are assumed to be the same in 
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(A3) and ( A4) , and the eventual results are independent of their 
numerical values. Except for an extremely thin laminar sub-layer, 
U is negligible compared to Ut . Then, in the turbulent region, 
(A3) and ( A4) substituted into (A2) gives: 


3 2 m _ 1_ ^m . 2 3 2 m 
3r 2 r 3r 3z 2 


where 


£ 2 = Ut i / Ut t = I 3m/ 3z I / I 3m/ 36 I 


( A5 ) 


As a first approximation we can assign a constant value to £. 
With £ constant equation (5) is solvable by separation of 
variables. The equation in z is an exponential: 


m( r, z) = m( r, 0) exp( z/ zo ) , where (A6) 

zo is the z decay constant whose numerical value comes 
from the solution of the equation in r. 

The equation in r is a Bessel equation with the following general 
solution: 

m( r , 0 ) = r^CaJi + C 4 Y i j > or 

m(r > 0) - = fca, ♦ c 2 Y, (Jf)l (A7) 

mu L v z d J 2 o J 

Cl = C 3 V.M 

C2 - c 4 V» n 

Jt is the Bessel function of the first kind, first order. 

Yi is the Bessel function of the second kind, first order. 

A. 2 BOUNDARY CONDITIONS 

The three unknowns Ci , C 2 , and zo require three boundary 

conditions. 

Boundary Condition (M. Since we have eliminated the laminar 

sub-layer from the solution, we will make the condition at the 

inside edge of that layer (r=r„) the first boundary condition. 

The "universal velocity profile" (3) places the edge of the 
laminar sub-layer at y* = 5, where 

y * = = 5 ( A8) 

Also in the laminar sub-layer, y* = v* where 


V 


v, (r „,z ) 


5 


( A9 ) 


Equation ( A9) introduces the wall shear, and the Blasius 
expression for wall shear [3] makes the connection back to the 
maximum angular momentum (which is the input quantity): 


T o = 0.0225 P v„ 2 f— ) *< Ren ) 

- Th 


-H 


Re» 


. Pv n r fl 


v« is the velocity at maximum angular momentum, 
r* is the radius at maximum angular momentum. 


( A1 0) 


Substitution of (10) into (9) gives Boundary Condi t i on : 

C l 'xl/ft — 1/ft 

— 1 j ( Re« ) (All) 

1 -rn J 

It should be noted that a guess of m must be used for the first 
trial solution, then succeeding trial values are obtained by 
setting the derivitive of the momentum equation ( eq. ( A1 3) ) equal 
to 0 and solving for r = r«. 

Boundary Condition 1t2 . Zero shear stress at the free surface. 

At the free surface turbulent fluctuations disappear; therefore, 

Ut • = 0. The condition, T r • = 0 at rt yields: 

( 8m/ 3r) . = 2m< r . , 0) / r . ( A1 2) 

l l l 

Equation ( A7) differentiated is: 

f 01 = — |"c3 Jo [~ ] + C4 Yo[~ ]] . CA13) 

3r 2 d L L za J L Zd J J 

Jo is the Bessel function of the first kind, 0 order, 

Yo is the Bessel function of the second kind, 0 order. 

Equations ( A1 2) and (A13) form the second boundary condition. 

Boundary Condition i\3. Maximum angular momentum. 

The maximum angular momentum, mn , is considered to be a given 
quantity. The angular momentum profile is normalized by dividing 
all m( r, z) by oih. 

Sol ut i on 


The solution for Cl, C2, and Zd (program ZGROT) depends upon the 
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fullness of the cylinder, r»/r 0 (ri/r 0 =0 being full, and ri/r«=1 

being empty) , and the Reynolds Number range, as outlined in the 

following table: 


TABLE A1 . For 10* < Re < 10® Re r . » - 5*10 7 



0. 01 

1 1 . 926 1 

3. 31 7E-7 

1 

0. 2748 

1 

0. 6609 

0. 2 

1 1.935 1 

4. 5 3 5 E - 2 

1 

0. 2730 

1 

0. 6628 

0. 4 

1 1 . 958 1 

0. 5564 

1 

0. 2584 

1 

0. 6823 

0. 6 

1 0. 9068 1 

2.175 

1 

0. 2065 

1 

0. 7363 

0. 8 

1 -0. 8738 1 

-3. 1 489 

1 

0. 1 229 

1 

0. 8374 

1 . 0 

1 1 

l 1 


1 

1 

0. 0 

1 

1 

1.0 


9. 323E + 8 
9. 336E + 8 
9. 477E+8 
9. 930E+8 
1 . 1 20E+9 


* Ft is a dimensionless quantity from which the wall shear can be 
calculated: 

Ft = Ti p( r«/U) 2 ( Re r i i / Re) 2 

Figure A1 shows the computed angular momentum profiles. 



angular momentum 



RADIUS RATIO, R/RO 


FIGURE A1. 


Program ZGROT also computes the factor S, the ratio of axial to 
radial angular momentum derivitives, and the average of Z over r 
once Cl, C2, and Zd/(roS) are computed. 
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TABLE A2 


r t /r 8 , 

0 . 01 
0 . 2 
0. 4 
0 . 6 
0 . 8 


Zi avq) 

1 . 42 
1 . 49 
1 . 53 
1 . 65 
1 . 70 


z d / r o 

0. 391 
0. 407 
0. 409 
0. 341 
0. 209 


These values of z d /r 0 indicate a large decay of angular momentum 
in the z direction. For example, in an axial distance of one 
radius the angular momentum or angular velocity would decay to 
1/10 its value at z = 0. A simple visual experiment shows this to 

be incorrect. 

A. 2 NUMERICAL SOLUTION OF THE MOMENTUM EQUATION WITH VARIABLE Z 
(Other boundary conditions and assumptions as in A. 1) 


In the linear solution of section A. 1 an average (constant) value 
of S was used. The radial variation of Z calculated from the 
resulting equation shows a spike in the region of m« , so that the 
assumption is inconsistent with the results. (This is 
paradoxical because a solution which allows a variable S results 
in a more nearly constant Z) . To test the possibility that the 
constant S assumption resulted in erroneously small values of 
2d/r fl) the following non-linear equation (which makes no 
assumptions concerning the momentum derivitives) was solved 
numerically. The equation contains the complete expressions for 
Tr • and T* z , equations ( 3A) and (4A) so that the computation 
should be valid into the laminar regions. ( To obtain the same z 
variation as before it is only necessary to assume similar 
prof i les. ) 

z variation (equation A6) . 


r variation: 


r 


3 2 m 
3r 2 




3 

3r 


3m 3m ^ m 
3r 3r z 2 


+ 2CR 2 - 
V 



0 ( A1 4) 


C and R are the "mixing length" parameters which were 
discussed and evaluated in the midterm report. 


A sample calculation for ri/r o =0.2 resulted in Zd/r o =0.464 which 
is still too small (observations would indicate 2d/r 0 > 2). 



APPENDIX B 


ADAPTATION OF THE SCHULTZ -GRUNOW INTEGRAL MOMENTUM METHOD 


The Schultz-Grunow case of interest is one in which the axial 
spacing between the disk and housing is considerably greater than 
the combined boundary layers. Contrary to the present problem, 
however, the length of cylinder between the rotor and housing is 
not considered enough to contribute to the drag, so it is 
necessary to extend the theory. It is also necessary to account 
for partially filled cylinders (ri/r»>0) and for the existence of 
a free surface as would occur in a low gravity environment. The 
presence of an actual spinning disk is immaterial since the flow 
in the liquid core does not depend on the means (such as a jet or 
spinning disk) used to induce the rotation. 

The approach taken is that the central core of liquid rotates as 
a solid body and the viscous effects and secondary flows take 
place in the thin boundary layers next to the solid and free 
surfaces. This description agrees basically with the visual 
pe r c e pt i on. 

B. 1 EXTENSION OF SCHULTZ-GRUNOW TO ACCOUNT FOR A LONG CYLINDRICA L. 
SECTION AND AN INTERIOR FREE SURFACE 

The boundary layer thickness 5 is taken to be that thickness of 
liquid which, being dragged along with the spinning disk, has 
exactly the centrifugal force to balance the drag due to the 
radial velocity component. The tangential drag in terms of 6 is 
that due to Blasius C 33 and the 1/7 power volocity distribution: 

To cos(e)=Cfp(r(Uo-Uf)) 7/4 (W6) lM ( B1 ) 

9 is the angle of the shear stress with respect to the 
tangential velocity. 

The next task is to determine 6, and that is accomplished by 
balancing the centrifugal force in the layer 6 against the drag 
force from ( B1 ) : 

f To si n< 9) 2Ttrdr = pruo 2 62nrdr ( B2) 

f = (contribution to radial drag force due to all of the 
surf aces) /( radi al drag force on the spinning disk 
alone) 

From ( B1 ) and ( B2) the boundary layer thickness is: 



6 = ( Cf tan( 0) ) 4 ' 5 r 3 ' 5 ( WUo ) ' ' 5 ( 1 -Ui /Uo ) 7 ' 5 . ( B3) 

The torque on one side of the disk is ( in our case integrating 
from rt instead of 0): 

r r ° 

To = r( 2ttr) To cos( 9) dr 
J r i 

Substitution from (B1) and ( B3) gives: 

To = C, puo 2 r, s ( 1 -U( / Uo ) 7 ' 5 0 -( r t / r, ) 2 3 ' s ] Reo ' 1 ' 5 

Where C. = 1.366 Cf 4/S (f tan( 6))' t/5 , 

Reo = Uo r s 2 / v. 

Equation ( B5) as developed is a more general case of the following 
Schultz-Grunow equation ( for one side of the disk) : 

To = 0.0311(1/2 puo z r 0 5 ) Re' w s ( B6) 


( B 4 ) 


( B5) 


C 9 can be evaluated from ( B5) and ( B6) by matching the conditions 
in ( B5) to those of ( B6) ; namely, ri=0, and Uf/Uo=0.5 (Schlichting 
[31), and f = 2 (end disk drag = spinning disk drag). Then, 

C, =0. 041 ( 2/ f ) 1 ' 5 < B7) 


rlindrical Section 


The tangential torque on one half of the cylinder is: 


Tc » i = 2nr „ c 2 LT« cos( 9) 

T. cos( 0) =C, p[ r. c Uf ) 7 ' * ( vt 6) 1 ' 4 

Assuming 6 the same as on the rotating disk we get: 

T f , i = 23 L rr.c] liM r ut/u. ] 7M 1 

To 5 r. [r, J Ll-Uf/UcJ 1 -( r i / r • ) 2 3 ' 5 

Fixed End Closure 

The torque due to one fixed end is: 
f r. 

Te = r( 2nr) T, cos( 9) dr , 

J r i 

and the shear stress is 

T. cos( 0) =Cf p( r U( ) 7 ' 4 ( v) 6) ‘ ' 4 . 

Again assuming the same 6 we get: 


( B8 ) 
( B9 ) 


( B 1 0) 


(B11) 
( B1 2) 
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It 

To 


( B1 3) 


[~ CJf/u)D 1 7 7 4 

Ll-Uf /U.J 

Inner Fixed Cylinder (optional) 

If an inner fixed cylinder exists, its torque contribution is the 
same as (BIO) with r,< replaced by r t : 

Xi _ 23 L |~r 1 1 1 * / 9 (Jf / (Jo .1 7 ' 4 1 (B14) 

To 5 r 0 Lr 0 J Ll-Uf/tdoJ 1-(ri/r 0 ) 

Since the torque on the disk (whether the disk is actual or 
imaginary) balances the fluid torque on all of the contacting 
surfaces, equation ( B5) gives the fluid torque. Changes in the 
rotation ratio u,/w D reflect the changes in geometry from that 
assumed in the original derivation. The following cases are two 
of the configurations of interest: 

Case 1 . For a partially or completely filled cylinder in zero 
gravity having a free surface ( ie. no torque) at ri> the 
following holds: 

1 = Tcyi/To + Te /To . 

This, combined with (BIO) and (B13) yields: 

wo / w, = L F( L) +1 J 4 ' 7 + 1 ( B1 5) 


where, 


F ( L) 


23 L r 

5 To L 


18/5 


1 - ( r » / r a ) 


2 3/5 


( B 1 6) 


Case 2 . Swirl in an open top vessel in a normal gravity 
environment: there will be no torque on the top end ( Te=0) . 

The following also allows for a possible inner fixed cylinder as 
in a proposed swirl test apparatus: 


1 = Tcri/To + Ti/To 


Wo / Uf 


[f( L) (l +( r i / r 0 ) 


t a / 3 


)] 


( B1 7) 


Now it is possible to calculate the factor f in equation ( B7) ( it 
should be pointed out that <2/f) ,/5 can be considered unity 
except for very long cylinders): 


f = 1 + f c + f e ( B1 8) 

If the liquid contacts a fixed end closure ft s 1; if not, as in 
case 2, fe=0. 



f 


1 

1-(ri/r fl ) l8/i 


( B1 9) 


c 


1 8 [ ~ Mf /(Jo 1 7 ' 4 L r 
5 Ll-Uf/UoJ r. L 



1 ■ / 3 


B. 2 EXPERIMENTAL CONFIRMATION OF THE ROTATION RATIO 

Since equations (B16) and (B17) involve observable rotation 
ratios, a visual experiment is useful for at least a crude 
check. For this purpose I constructed a plexiglass cylinder in 
which swirl was induced by a vertical axis spinning disk. 
Near-neutral-buoyancy particles made the liquid rotation 
visible. Both the liquid and disk rotation rates were timed by 
stopwatch. The specifics are: 


Fluid water 

ri 0 

r. 2. 27 in. -- 0. 0576 m 

r, t 2.13 in. -- 0.0544 m 

L varied from 0. 5 to 5. 625 in. 

Uo 3.5 rev/s (approx, constant) 


Figure B1 shows the results of 30 measurements at 6 values of 
L/ r. . 



Figure B1. Disk/liquid rotation frequency ratio as a function 
of cylinder length. 

The vertical bars show the approximate spread of the data. The 
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main experimental difficulty was in following particles which 
tended to drift in and out of the main rotating core. As a 
result, the timing of rotation was imprecise. 

The solid line represents the theory. From equation (B16) 

F( L) = 5. 76( L/ r. ) , 

and the rotation ratio Ud/u* was obtained from equation ( B1 5) . 


These simple experiments cannot actually check the theory since 
the frictional torque was not measured; however, the effect of 
the added cylinder length can be checked since length is 
reflected in the change in u«/Ud from the Schultz-Grunow case 
where / Uo =0. 5. Even though the observed Uo/U{ appear to be 
about 9 % higher than the predictions it seems reasonable to 
conclude that the observations generally substantiate the 
modifications. The frictional drag has already been compared 
with experiments as documented in Schlichting [ 3] , and the 
Schultz-Grunow predictions were found to be 17% low over a 
Reynolds Number range 2(10) 5 to (10) 7 . 
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